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We use the Hellman-Feynman (HF) and Hypervirial (HV) theorems, to 
calculate the perturbative coefficients of the eigenenergies formal series, 
in the case of the Coulomb potential with a radial linear term and the 
f> i 1 radial quartic anharmonic oscillator potential. This calculation method, 

contrary to the usual Rayleigh-Schrodinger Perturbation Theory (RSPT), 
does not require the calculation of eigenfunctions coefficients. This method 
is a fast and efficient tool for the calculation of large order eigenenergies 
coefficients. 
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Abstract 



1 Introduction 



Many physical phenomena encountered in different fields of physics, require the 
modelization of the interaction through potentials that are believed to account 
for the physical mechanisms occurring in the system. The construction of these 
potentials is constrained by physical consideration in their analytical forms and 
introduces parameters. For instance, the problem of the electron in the field 
of a nucleus is tackled with a Coulomb potential, since it is known that this is 
the interaction which is in work. In nuclear physics, the nuclear energy levels 
are obtained by modelizing the interaction, the nucleon feels through (e.g.) a 
Woods-Saxon potential. In particle physics, the mesons are described as quark 
antiquark bound states through a phenomenological model where the potential 
is given at once. 

While bound states problem can be solved analytically [1] for few potentials, 
among which Coulomb and harmonic oscillator potentials, for most of poten- 
tials of physical interest, one has to resort to direct numerical techniques or to 
approximation methods. Note that, beyond the physical problem itself, these 
solving techniques have become a field of study in their own. 
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Perturbation theory is one of the most widely used approximation method to 
solve eigenvalue problems in theoretical physics [2, 3, 4]. Mathematical aspects 
of perturbation theory are treated in monographs [5, 6, 7], while applications in 
various fields are described in books [8, 9, 10]. 

Rayleigh-Schrodingcr Perturbation Theory (RSPT) is designed for systems 
whose Hamiltonian can be partitioned into a completely solvable part with 
known eigenvalues and eigenfunctions, and a perturbation part [2, 3, 4, 11]. The 
formalism of RSPT expresses the eigenenergies and the eigenfunctions as formal 
series in the choosen perturbation parameter. Accurate results for physical ob- 
servables have been obtained from this perturbative method. However, one of 
the shortcomings of the RSPT approach is that it requires, in general, consider- 
able computational time and effort. Indeed, this method uses the eigenfunctions 
solutions which are time consuming and intermediary results in an eigenvalue 
problem [4, 11]. 

For various interactions, a way to bypass the computation of the eigen- 
functions is to use Hcllman-Feynman (HF) and the Hypcrvirial (HV) theorems 
[12, 13] to generate eigenenergies perturbative coefficients [14, 15, 16, 17, 18]. 

The purpose of the paper is to use these theorems (denoted by HFHV 
method) to calculate the eigenvalues perturbative coefficients for two poten- 
tials of physical interest namely: the Coulomb potential plus radial linear term 
and the radial quartic anharmonic oscillator potential. The eigenenergies coeffi- 
cients are computed recursively without the need to calculate the eigenfunctions. 
Moreover, we show that the application of the HFHV method requires consid- 
erable less computational time comparatively to the usual RSPT method and 
allows a fast and efficient calculation of very large order eigenvalues coefficients. 



2 Method of calculation 

2.1 Hellman-Feynman and Hypervirial theorems 

Let us recall Hellman-Feynman (HF) and Hypervirial (HV) theorems. The HF 
theorem states that if we consider a system with a Hamiltonian H = H(X), 
depending on some parameters A, with normalized eigenvector \f(A) and eigen- 
value E(X), then: 



!£-<•<*> 



OH 



*(A)> = <^> (1) 



Let H be a timc-independant Hamiltonian with a normalized bound state 
eigenfunction and Q a time-independant linear operator. If the expectation 
value of 6 over the stationnary eigenstate * is not infinite, then the HV theorem 
states that, for the operator 6: 

(*|[e,fl]|*) = <[e,fl]) = o (2) 

where [0,-ff] = OH — H is the commutator of H and O. Physically, this, 
is just the statement that, for a stationary state, the expectation value of is 
time-independant. 
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2.2 General formalism 



The two-body stationnary Schrodinger equation with central potential V(r), 
can be written (in h = c = 1 units ) as: 



1 d 2 . , 1(1 + 1) 

h V(r) + — 

2fj,dr 2 K > 2(j,r 2 



Vn,l(r) = (3) 



where /i is the reduced mass of the system, E n j and y n ,i(f) are respectively the 
eigenvalue and the reduced radial wave function, describing a bound state of 
radial and orbital angular momentum quantum numbers n and / respectively. 
The wave function obeys the following boundary conditions: 

VnM = l/n,l(°°) = ( 4 ) 

Setting x — Xr, the scaling property of the Schrodinger equation allows to write 
(3) as: 

H y n ,l(x) = e n .i y n ,l( x ) ( 5 ) 

with: 

H = -^ + U( X ) (6) 

and: 

U(x) = V (x) + + g V 1 (x) (7) 

x z 

The perturbation is cast in the term g V\(x) where g is the controlling parameter, 
i.e., the expansion parameter. 

In the usual RSPT method, the eigenvalues and eigenfunctions are written 
as series expansions in the perturbative parameter g: 

oo 

e n ,i = J2 e ^(k) 9 k (8) 

k=0 

oo 

Vn,l( X ) = ^2Vn,l( k ' X ) 9 k (9) 
fc=0 

where c n ,l(k) and y n ,i(k,x) are, respectively, the eigcncnergies and eigenfunc- 
tions perturbative coefficients, k — in the expansion corresponds to the unper- 
turbed system, i.e., e ni ;(0) and y n j{0,x) are the eigenvalues and eigenfunctions 
for the potential (7) with g = 0. 

With the RSPT method, the calculation of e n ,i{k) coefficients requires then 
the knowledge of y n ,i(i,x) where i — 0, • • • k — 1 [2, 3, 4, 11]. This is the 
shortcomings of the method. 

It happens that, for various interactions, and thanks to HF and HV theorems, 
the e„ ; ;(fc) coefficients can be given by a recursion formula (see below) and thus 
there is no need to calculate the intermediary eigenfunctions [14, 15, 16, 17, 18]. 

2.3 The potentials 

We consider two central potentials of physical interest, the Coulomb potential 
with a radial linear term and the quartic anharmonic oscillator potential. 
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The Coulomb potential with a radial linear term: 



V(r) = --+0r a,(3>0 (10) 
r 

corresponds to a spherical Stark effect in hydrogenic atoms and is, also, a suc- 
cessful phenomenological model for the description of heavy quark antiquark 
bound states in particle physics [19, 20]. In this case, taking A = 2/ia, the 
eigenvalues are given by: 

a 2 2n 

and the scaled potential (7) reads: 



e.J = 4=^ (11) 



where g is the dimcnsionessless parameter used in the expansion. The term 
gx is thus treated as a perturbation, while the the unperturbed Hamiltonian 
is the usual radial Coulomb Hamiltonian whose eigenfunctions and eigenvalues 
are known [1, 11]. 

The second potential is the quartic anharmonic oscillator potential: 

V(r)=ar 2 +br 4 a, b > (13) 

which has received a great deal of interest because of its importance in the un- 
derstanding of molecular vibrations [21], in some areas of solid state physics 
[22, 23] and in Quantum Field Theory [24, 25]. Furthermore, the quartic anhar- 
monic oscillator serves as a testing tool to check various approximative methods 
in Quantum Mechanics [26, 27]. Setting A = (2/x a) 1 / 4 , the eigenvalues are given 
by: 

eml = \—E n j (14) 
V a 

and the scaled form of this potential (7) reads: 

where g is the perturbative dimensionessless parameter. Here again, the term 
gx 4 is treated as a perturbation, while the unperturbed Hamiltonian is the 
radial harmonic oscillator Hamiltonian whose eigenfunctions and eigenvalues 
are known [1, 11]. 

2.4 Recurrence equations on the e n:l (k) coefficients 

Let us consider the time-independant linear operator as 6 = x^ 1 -^, where 
j is an integer such that j > [14, 15, 18, 28, 29]. Taking the Hamiltonian in 
the form (6), one has: 
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From the definition (5), one can easily show that 1 : 



(at ^) = ( x* U(x) ) - e n ,, (a? ) (17) 



Integration by parts, using the boundary conditions (4) together with the asymp- 
totic behavior y n j(x) <~ x l+1 for small x, leads to: 

{x^ 1 ^) = -\{j-l){^- 2 ) (18) 

Inserting relations (16-18) into equation (2), the HV theorem provides relation- 
ships between e„ ; ; and the various expectation values (x^ ) through the following 
equations [14, 15, 18, 28, 29]: 

2(3 + 1) e nil (X* ) = -\j(f-l)( x>- 2 )+2(j + l) ( x? U(x) ) + ( x* +1 ) 

X (19) 

The HF theorem (1) allows to derive a further relation: 



de n ,i = , dU(x) 
dg dg 



\ - l \r ( \ \ - J ( x ) ior the potential (12) 1 
( x i) for the potent i a i (15) J 



The essence of the application of the HF and HV theorems (HFHV perturbative 
method) is to assume, in the spirit of perturbative methods, that the energy 
and the expectation values of position coordinates can be expanded in power 
series of the perturbation parameter g as (see eq. (8) for the energy): 

00 

<^>=5>f g k (21) 

fc=0 

By equating like powers of g on both sides of equations (19-20) after substitution 
of expansions (8) and (21), we find a set of coupled relations involving e nj i(fc) 

and x^ . 

To be more precise, let us consider the Cornell potential defined by (10) or 
(12). In this case, equation (19) can be written as: 

e n>1 (x> )= aj ( x>- 2 xi- 1 ) + g 7, ( ) (22) 

with: 



((2t + l) a -j a ) . 2j + l. 2_±±l m) 

a '~ J 4j + 4 ' 2j + 2' 7 ^ 2. ? + 2 

Inserting the expansions (8) and (21), in equations (22-20), the identification of 
the coefficients of powers of g leads to (j integer > 0) : 



J2^i(k-P) =a j xf} 2 +(3 j xf} l +e{k-l) lj xf- l 1) k>0 (24) 

(fc-i) 

e n ,l(k) = ^7— * > 1 (25) 



1 Throughout this paper, the notation < £ > stands, the quantum number indices n and 
I being understood, for the expectation value < £ >= (\P | £ | Nt), where \P represents the 
normalized eigenstate of the hamiltonian H under consideration. 
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where 8(t) is the Heaviside Step function. It is clear from the above equations 
that the starting point is the normalization condition: 

4 0) = 1 and 4* =° for fc > 1 (26) 
and the Coulomb eigenvalues: 

^ = - 4(n+i + l)» (2?) 

The set of equations (24-25) together with (26-27) allows one (see below) to 
calculate all the e„ ; ;(fc) coefficients. 

In a similar way for the quartic anharmonic oscillator (13) or (15), taking j 
even (j = 2i), equation (19) reads: 

e n ,l (x 21 > = a, < x 21 ' 2 ) + ft ( x 2 *+ 2 )+ gil ( x 2l+i ) (28) 

with: 

, ((21 + l) 2 - 4i 2 ) 2i + 2 2i + 'i 

4TT2 5 —r 7<= 2TTT (29) 

By analogy with equations (24-25) and noting x^ = x\ p \ one obtains in the 
quartic anharmonic potential case: 

k 

£>„,,(*: -p) x\ p) =a i x\ k J 1 +p i x\%+0(k-l) 7i x? + - 1) fc>0 (30) 

p=0 

~(fc-l) 

e n Ak) = ^ V - k>l (31) 

together with the normalization condition: 

4 0) = 1 and x { k) = for fc > 1 (32) 
and the harmonic oscillator eigenvalues: 

e„,i(0) = 3 + 4n + 2Z (33) 

2.5 Calculation of the e n> j(fc) coefficients 

To calculate explicitely the e n ,i{k) coefficients from the equations (24-25) or 
(30-31), one proceeds in a hierarchical manner by giving to k various integer 
values starting with fc = 0. Let us describe this method in details for the case 
of the Cornell potential. Consider equation (24) for fc = 0, one obtains: 

e n ,,(0) xf =a j xf} 2 +f3 j xf} 1 (34) 

which is a recurrence of depth three on the coefficients. One needs two 

initial conditions to calculate all the x^p coefficients. One condition is given 

by the normalization (26), namely Xq ^ = 1, the secund one is given by taking 
j = in cq. (34) and one has: 

- - 2 e -'<°> - wtttw (35) 
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The recurrence equation (34) together with the initial conditions (26) and (35) 
allow to calculate all x^ coefficients. The next step is to consider equation (24) 
for k > 1, which can be written, using (25), in the following form: 

en,i(0) xf } =a j xf} 2 +p j xf} 1 +S(j,k) k>l (36) 

where S(j, k) depends only on x^ terms with p < k and reads: 

fe-i ~(fe-i-p) Jp) 

5(.%k) = -J2 1 . *' + 7^ 07) 
p=0 ^ 

The obtained equation (36) is an inhomogenous recurrence of depth three on 
the Xj k ^ coefficients (k kept fixed), S(j,k) being the inhomogenous term. In 

(k) 

the same way, one needs two initial conditions to calculate the Xj coefficients. 

One condition is given by the normalization (26), namely x^ — (k > 1), the 
secund one is given by taking j — in eq. (36) and one has, using (26): 

x « = ( 3 -£) 4 fe - 1} k>l (38) 

One proceeds step by step by giving to k various integer values in equations 
(36-38), starting with k = 1. 

Finally , we can, thus, calculate the r^ 1 pcrturbative coefficient e n j(r ) 
through the relation (25) from the knowledge of the xf 1 with 2 < k < ro — 1 
and — 1 < j < ro — k. 

For the quartic anharmonic oscillator (15), similarly, one obtains (fc = 0): 

^^=6^(0) xf-aafl^ (39) 
4° } = 1 and xf ] = e„,j(0)/2 = 3/2 + 2n + I (40) 

Pi 4% = en,i(0) ^-ai^+^i,*) (41) 

fc-l -(fe-l-p) Mp) 
8{i, k) = - 7 < x^ + ^ 2 . *' (42) 

*f = and if = (-3 + yfc) s<*-i) (43 ) 

Consequently, the rg h perturbative coefficient e„.;(ro) is calculated through the 
relation (31) from the knowledge of the xf with 3 < k < r — 1 and < j < 
ro-k + 1. 



and (fc > 1): 



2 For j = —1,0, the x^' arc given by the initial conditions (26, 35 and 38) and the recurcnccs 

Ik) 

displayed above allow to calculate the Xj ' for j > 1. 

3 For j = 0, 1, the x^ arc given by the initial conditions (40 and 43) and the recurcnccs 

^(k) 

displayed above allow to calculate the x\ for j > 2. 
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3 Results and discussion 



We have implemented the HFHV perturbative method described above for the 
calculations of e n ,i{k) coefficients, at any order k, into a Mathematica routine. 

The first eight eigenergies coefficients are displayed in Tables 1-2 for both 
the considered potentials. 

Tablel. The coefficients e n j(k) in the case of the Cornell potential for few 
states. 



State 


n = 0,Z = 0(15) 


n=l,l = 0(25) 


n = 0,1 = 1(1P) 


k = 


i 

4 


i 

1« 


i 

in 


k = 1 


3 


12 


10 


k = 2 


-12 


-528 


-480 


k = 3 


216 


105984 


99840 


k = 4 


-6360 


-34775040, 


-33054720 


k = 5 


245952 


14770372608 


14051573760 


k = 6 


-11433984 


-7410570756096 


-7038283284480 


k = 7 


610773696 


4197622600433664 


3976876146032640 



Table2. The coefficients 4 fe e n ,i(k) in the case of the quartic anharmonic 
oscillator for few states. 



State 


n = 0,l = 0(15) 


rc = 1,1 = 0(25) 


n = 0,1 = 1(1P) 


k = 


3 


7 


5 


k = 1 


15 


75 


35 


k = 2 


-165 


-1575 


-525 


fc = 3 


3915 


66825 


16625 


k = 4 


-520485/4 


-15184575/4 


-2894325/4 


fc = 5 


21304485/4 


1024977375/4 


152440575/4 


k = (5 


-2026946145/8 


-155898295875/8 


-18353729625/8 


k = 7 


108603230895/8 


12977225578125/8 


1224596281125/8 



Table3. CPU time (in secunds) for the calculation of eo,o(&) (k — 0, ro) 
coefficients in the case of the Cornell potential for the IS state. 





ro = 10 


ro = 20 


r = 30 


r = 40 


RSPT method 


26.2 


188.7 


1024.9 


2626.2 


HFHV method 


0.016 


0.032 


0.062 


0.16 



Table4. CPU time (in secunds) for the calculation of eo,o(k) (k = 0, ro) 
coefficients in the case of the quartic anharmonic oscillator for the IS state. 





ro = 10 


r = 20 


r = 30 


r = 40 


RSPT method 


43.4 


277.8 


1213.6 


2676.6 


HFHV method 


0.016 


0.047 


0.11 


0.27 



We have also compared the present method to the usual RSPT one 4 , in the 
simple case of the IS state. For both potentials, the respective CPU times 5 are 
displayed in Tables 3-4. One can sec from the latter tables that the computa- 
tional time of the HFHV method is reduced by a factor 10 4 comparatively to the 

4 For the RSPT method, we have implemented, in a Mathematica program, the set of 
equations displayed in section 7.3 of Ref. [4]. 

""Throughout this paper, the mcntionncd CPU times correspond to runs on a PIV personnal 
computer with 1.5 Go of RAM memory. 
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usual RSPT method. In addition, let us note that the ratio of the respective 
CPU time (RSPT/HFHV) grows with increasing order. 

Furthermore, the HFHV method reveals itself as a fast and efficient tool for 
a systematic calculation of large order perturbative coefficients. Indeed, consid- 
ering, for example, the IS state in the Cornell potential case, the necessary CPU 
time to calculate the first 1000 perturbative coefficients is less than four (04) 
hours 6 . Consequently, one can find, in principle, the asymptotic behavior of the 
£n,l(k) coefficients for large k, which reads, for both the considered potentials: 

e n ,l( k ) ~ (-l) fe+1 T(k + b(n,l)) a(n,l) k for k -» oo (44) 

where the notation T stands for the Euler Gamma function. For example, for the 
IS state, considering the Cornell (resp. quartic anharmonic) interaction, we have 
found that a(0,0) ~ 6 (resp. a(0,0) ~ 3/2 ) and 6(0,0) ~ 2 (resp. 6(0,0) ~ 3/2 
) 7 . Let us note that, for the quartic anharmonic potential, the asymptotic 
behavior (44) together with the numerical values of a(0, 0) and 6(0,0) are in 
perfect agreement 8 with the one dimensional calculation of Bender and Wu in 
the framework of WKB analysis and difference-equation method [24, 30, 31]. In 
the Cornell interaction case, the large k behavior (44) agrees fairly well with 
the WKB result [32]. 

Let us make some comments about the asymptotic behavior (44) of the 
e n ,i(k) coefficients. This large k behavior means that, for the two potentials un- 
der consideration, the expansion (8) has its radius of convergence equal to zero 
and is, thus, a divergent serie for all g. The divergence of the perturbative serie 
(8) indicates that taking, for the Cornell (resp. quartic anharmonic) potential 
the term gx (resp. gx 4 ) as a perturbation leads to a singular perturbation. For 
the Cornell (resp. quartic anharmonic) interaction, a simple way to understand 
this singular behavior is to compare, exp(— 2( - K _j^ +1 - ) ) (resp. exp(— \)), the 
controlling factor of the large x behavior 9 of the unpertrurbed (g — 0) eigen- 

function with, cxp(— 2g ■ ' x ' ) (resp. cxp(— 9 ' ^ x )), the controlling factor of 
the large x behavior for g ^ 0. One can see an abrupt change in the nature of 
the eigenfunction when we go, for large x, to the limit g — > 0+. For the Cornell 
(resp. quartic anharmonic) potential, this feature occurs because, for large x, 
the perturbation term gx (resp. gx 4 ) is not negligeable comparatively to —1/x 
(resp. x 2 ). 

The other remark deals with the following point. Since the the expansion (8) 
is a divergent serie, one may worry about the interest of the c n ,l(k) coefficients 
calculation. In fact, different mathematical methods (like the Pade approxi- 
mants) allow to sum divergent series and have been used with some success in 
different fields of physics [4, 33, 34, 35, 36, 37, 38, 39, 40, 41]. For the eigenergies 
expansion (8) in the case of the considered potentials, the choice of an optimized 
summation method is beyond the scope the present work. Anyhow, the interest 

6 For the other states, the CPU time has the same order of magnitude. 
7 These results have been obtained with a satisfactory accuracy from the first 1000 pertur- 
bative coefficients. For excited states, one needs more perturbative coefficients. 

8 Up to rescaling of the quantities g and e n> ; to match the definitions given by Bender and 
Wu. 

9 The notion of controlling factor for large x comes from the fact that x = oo is an irregular 
singular point of the dimcnsioncssless Schrodinger equation (5) with potentials (12) and (15). 
For more details, see Ref. [4]. 
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of the t n ,i{k) coefficients calculation lies in the fact that the accuracy of the ap- 
proximative sum increases with increasing known terms in the expansion, this 
feature being independent of the particular choice of the summation method 
used. 

4 Comments and conclusion 

Before concluding, some comments are in order. We have shown the uselfulness 
of HFHV method for Coulomb plus linear term and the quartic anharmonic os- 
cillator potentials that interest us for other issues. This method can be applied, 
mutadis mutandis, to radial interactions of the form V(r) = — 1/r + g r p (p > 2) 
or V(r) = r 2 +g r 2 ? (p > 3). 

More generally, the HFHV method has been used, in fact, for many interac- 
tions, where the potentials take the form of a power serie (or a finite polynomial) , 
the unperturbed Hamiltonian being the Coulomb one or the harmonic oscillator 
one (for a general discussion see Ref. [29]). For instance, the method was ap- 
plied for: the screened Coulomb potential [15], the TV— dimensional generalized 
exponential cosine-screened Coulomb potential 10 [42, 43], the potential V(r) = 
r 2 + Ar 2 /(1 + g r 2 ) [44], the Gaussian potential [45, 46], etc. More recently 
[18], it has been applied to a set of potential like the Gaussian potential, Patil 
potential, etc. To be complete, let us note that the HFHV method has been 
also extended to the relativistic case for the Klein-Gordon equation, with a 
Coulomb potential plus polynomial perturbation, to calculate the eigenvalues 
perturbative coefficients [17]. 

The interesting and powerful side of the method stems from the fact that one 
does not need the eigenfunctions. Indeed, expressing the potential as a power 
serie (or a finite polynomial) in the coordinates x, the HF and HV theorems 
provide recursion relations that give the eigenenergy coefficients through the 
knowledge of the expansion coefficients of the moments (x 3 ). The coefficients 
x^ are computed by column, in a hierarchical manner, beginning with xf \ 

In conclusion, by applying the HF and HV theorems to the interactions (10) 
and (13), we have shown that the eigenvalues perturbative coefficients e„ ; ;(fc) 
can be obtained, at any order k, without any wave functions calculations. Thus, 
the HFHV method provides a remarkably fast and reliable tool to calculate very 
large order eigenenergies perturbative coefficients. 

Acknowledgement We would like to thank S. Hassani for discussions and 
critical reading of the manuscript. 
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